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Abstract. The problem of Finding Achievable Region (FAR) is a mean- 
ingful and novel problem, which is overlooked by previous researchers. 
Informally, given there are a set of obstacles in 2D space and a maximum 
path length (MPL), one (e.g., a person) walks freely from a starting point 
and stop to continue walk once the total travelled-distance is equal to 
MPL, we call the location (s)he stops as terminate point, the FAR prob- 
lem is to find the union of all possible terminate points. In this paper, we 
first propose and formulate this problem, followed by presenting a slice 
based method for tackling this problem. We prove the correctness of our 
proposed method, and demonstrate that the method takes 0(A 3 log A) 
time for the FAR problem. 



1 Introduction 

The problem of Finding Achievable Region(FAR), viewing from different per- 
spective, is similar to the following several classical problems, Visibility Polygon 
problem pQ, Visibility Graph problem [2J, Art Gallery problem [3], Watchman 
Route problem [3], Sensor Coverage problem |5I6) and the Motion Planning 
problem [7J, etc. For these classical problems, which have been well studied and 
a large bulk of solutions can be found in recent publications |8l5j . The FAR 
problem is a novel problem, which is overlooked by previous researchers. For 
ease of discussion, we, in here, present a visual example to illustrate what the 
FAR problem is. 

The black point is the starting point, the maximum path length (MPL) is 
the radium of the circle, and we assume the black line segment is an obstacle 



as shown in Figure 1(a) Figure 1(b) illustrates the result of FAR, i.e., all the 
locations in the circle are achievable except the locations in the white region. 
Therefore, the achievable region (AR) is the subtraction between the circle and 
the white region. Note that, in this example, there is only one obstacle, in this 
case we can obtain the result easily. However, once there are a large number of 



obstacles like Figure 1(c) illustrates or more complicate cases, the FAR problem 
will turn into rather complicate. In this paper, our goal is to develop an accurate 
and efficient method for tackling the FAR problem. 

The FAR problem can be took as a promising problem since it can be mapped 
to many application scenarios. For instance, assuming that a criminal gangs 
robbed a bank at 3 a.m. and now they are running away by riding "Toyota 



Royal V6", then, based on the geographical information around the bank, the 
maximum speed of the criminal car and the elapsed time, the policeman can 
get the achievable region the criminals may locate, which will contribute to the 
reduction of search range, note that, a more tight search scope will save more 
search cost in terms of manpower and time; or assuming that there is a lost 
child, someone told his mother that (s)he have seen the child in position "A" 
an hour ago, then, based on geographic information around location "A" and 
the maximum speed that the child can walk, the mother can get the achievable 
region her child may locate, then she may turn to her relatives for help and 
seek the lost child in the achievable region. In our real life, there are many 
application scenarios like the above two examples, we do not further enumerate 
them for saving space. People possible argue that the criminal car or the lost 
child is continuous moving, and the obstacles possible are various geometries 
in addition to line segments. Certainly, it is needed to consider these factors 
in real applications. It is worth noting, the examples we presented above just 
illustrate the potential application scenarios for the FAR problem. In addition, 
line segment is the most elemental element for various geometries. In this paper, 
we are interesting in line segment obstacles. 




(a) initial stage (b) achievable region (c) complicate case 

Fig. 1. Example of Finding Achievable Region 



For tackling the FAR problem, there are two main challenges. Data represen- 
tation for achievable region as well as intermediate result is a vital issue that we 



should make well consideration. Taking Figure 1(b) for example, the achievable 
region is the subtraction between the big circle and the white region, then, how 
to represent the result of this subtraction, or how to represent the white region? 
Another main challenge is the design of algorithm, that is, how to develop an 
algorithm that can efficiently and accurately find out the achievable region? 

Motivate by the above presentation, we, in this paper, investigate the FAR 
problem. To sum up, we make the following contributions. 

— We first propose and formulate the FAR problem. 

— We delve into the details of this problem, and develop an efficient and accu- 
rate approach for tackling it. 

— We, through theoretical analysis, prove the correctness of our proposed 
method and demonstrate its' computational complexity. 
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The rest of the paper is organized as follows. In Section [51 we review related 
works. We, in Section [21 formally formulate this problem and introduce the 
main notations used in this paper. In section [4l we present a slice based method 
for solving the FAR problem. In particular, we prove the correctness of our 
proposed method and analyse its computational complexity. Finally, we conclude 
this paper and present our future research directions in Section [5] 

2 Related Works 

There are several classical problems (such as visibility polygon problem and 
Sensor coverage problem) which are most related to the FAR problem. To some 
extent, these problems can be classified into two categories: computation of vis- 
ibility and coverage problem. 

2.1 Computation of Visibility 

Visibility is a mathematical abstraction of the real-life notion of visibility. Given 
a set of obstacles in the Euclidean space, two points in the space are said to be vis- 
ible to each other, if the line segment that joins the two points does not intersect 
any obstacles. Computation of visibility is among the basic problems in com- 
putational geometry and has many applications in computer graphics, motion 
planning, and other areas [5]. Specifically, there are several classical problems 
which can be view as the variants of visibility problem. 

One of variant is the visibility polygon problem, in particular, the visibility 
polygon or visibility region for a point, say p, in the plane among non-transparent 
obstacles possible is an unbounded polygonal region, which is consisted of all 
points of the plane visible from p. On the other hand, if the visibility polygon 
is bounded then it is a star-shaped polygon ilj. The visibility graph problem 
is another variant, a visibility graph is a graph of inter-visible locations, typi- 
cally for a set of points and obstacles in the Euclidean plane. Each node in the 
graph represents a point location, and each edge represents a visible connection 
between them [2]. Moreover, the art gallery problem or museum problem also 
is a well-studied visibility problem, which originates from a real-world problem 
of guarding an art gallery with the minimum number of guards which together 
can observe the whole gallery [3]. The watchman route problem is an optimiza- 
tion problem where the objective is to compute the shortest route a watchman 
should take to guard an entire area with obstacles given only a map of the area. 
The challenge is to make sure the watchman peeks behind every corner and 
to determine the best order in which corners should be visited in. There are 
polynomial-time solutions but they all suffer from severe numerical problems in- 
herent in the computations [4 . Last but not least, the motion planning problem 
is to produce a continuous motion that connects a start configuration S and a 
goal configuration G, while avoiding collision with known obstacles. The robot 
and obstacle geometry is described in a 2D or 3D workspace, while the motion 
is represented as a path in configuration space [7]. 
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Beside the above problems can be classified into the computation of visibility, 
fortress problem [5] , hidden line removal problem [TO] and hidden surface removal 
problem [11] can also be classified into this category. Obviously, all these problem 
are different from the FAR problem. 

2.2 Coverage Problem 

The term "coverage" is used in many fields. For instance, fault coverage in VLSI 
[12 , fold coverage in genetic sequencing [T3], code coverage in software testing 
[14] . etc. In here, we mainly refer to the coverage problem in sensor networks. 

Coverage which is one of the most important performance metrics for sensor 
networks reflects how well a sensor field is monitored. Individual sensor coverage 
models are dependent on the sensing functions of different types of sensors, 
while network-wide sensing coverage is a collective performance measure for 
geographically distributed sensor nodes |5l6j . Though coverage problem in sensor 
network and the FAR problem share some similarities since both of the problems 
involve the concept of "region" , it is not difficult to find that they are different 
in essence. 

3 Problem Definition and Notations 

At frist we give the definition of path length before we formally define the problem 
of finding achievable region. 

Definition 1. Given a starting point, SP, and a terminate point, TP, in 2D 
space, one freely walks from SP to TP, the path length is the total travelled- 
distance from SP to TP. 

Note that, if the SP and TP is the same point in 2D space, the path length 
is not definitely equal to 0. For example, given two different locations, say A 
and B, in 2D space, and we assume that both SP and TP are in location A, if 
one person firstly walks from A to B, and walks back from B to A, then, the 
length path is larger or equal to 2 * ||rfis(j4, B)\\, where ||<iis(yl, B)\\ denotes the 
Euclidean distance between the two locations A and B. 

Definition 2. Given there are a set N of obstacles in 2D space. For any two 
obstacles Oi and Oj (i =fi j, 1 < i, j < N), we assume Oi \\ Oj (i.e., the intersec- 
tion of the two obstacles is %). Given a starting point, SP (SP £ |J^. 1 Oi), the 
maximum value of path length is MPL. One starts (from the SP) to freely walk 
but cannot directly pass through any obstacle, and stop to continue walk once 
the travelled-distance is equal to MPL. The problem of Finding Achievable 
Region is to find all the possible locations (or points) it can reach at when SP 
and MPL is given beforehand. 

Note that, the obstacles, in real-world, possible are various geometries such as 
rectangle, circle, ellipse, trapezoid, line segment, even to other irregular shapes. 
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Table 1. Main symbols used in this paper 



symbols 


Description 


N 


number of obstacles 


Oi 


the i obstacle 


CAR 


candidate achievable region 


COs 


candidate obstacles 


AR 


achievable region 


MPL 


maximum path length 


Sj 


the j slice 


AARj 


accurate achievable region in Sj 



Since line segment is the most elemental element for other geometries, we, in this 
paper, only consider the obstacles are line segments, it will lay a foundation for 
further investigating the generic case, that is, the obstacles are arbitrary shapes. 
In the rest of this paper, if without stating otherwise, we use obstacle(s) to 
substitute line segment obstacle(s), and we assume we have indexed all obstacles 
in a data structure, such as R-tree, based on their MBRs. For convenience, we 
summarize the main symbols in Table [T] 

4 Slice Based Method 

For tackling the FAR problem, a slice based method is proposed, we first intro- 
duce the framework followed by presenting the details. 



4.1 Framework 

The framework of our solution is illustrated in Figure [5] Specifically, we first let 
the circle centering at SP and with radium MPL as candidate achievable region 
(CAR), and use the CAR to prune those unrelated obstacles and unrelated part 
of some obstacles, thus obtain the candidate obstacles (COs) (line 1-2). Next, 
we divide the CAR into multiple parts based on slicing rules (line 3). Then, we 
evaluate accurate achievable region ( A AR) for each slice Sj , after we dealt with 
each slice, we combine all the AARs from each slice and return the result at last 
(line 4-7). 



4.2 Pruning Unrelated Obstacles 

It is easy to know that AR is equal to CAR when there is no any obstacle, and 
AR must be a subset of CAR (i.e., ARC CAR) once there exist obstacle(s). 
For instance, in Figure |T(b)j the non-white region (corresponding to AR) is a 
subset of the big circle (corresponding to CAR). Therefore, at the beginning 
of evaluating the AR, we can use the CAR to prune some unrelated obstacles. 
Specifically, there are three types of obstacles that should be pruned or be pruned 
partially. 
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Procedure Framework.FAR { 

Input: J2iLi Oi, M PL, SP 
Output: AR 

(1) AR <- 0, CAR. «— circle centering at SP and with radium MPL 

(2) prune unrelated obstacles and part of obstacles 

(3) divide CAR, based on slicing rules, into k slices, Si, S2, • • • , Sk 

(4) For j = 1 to k 

(5) AARj evaluate the accurate achievable region in Sj 

(6) AT?, <- U 3 fe =1 AAA, 

(7) return AR } 

Fig. 2. Framework of Finding Achievable Region 



Lemma 1. Given there are a number N of obstacles in 2D plane, let circle 
centering at SP with radium MPL be the CAR. Any obstacle Oi (1 < % < N) 
can be safely pruned if Oi ft CAR ( i.e., Oi and CAR are disjointed). 

Proof. We only need to prove any sub-segment of the obstacle Oi makes no im- 
pact on the AR under the assumption that Oi \] CAR. Without loss of generality, 
we cut Oi into m segments, say si, S2, ••• , s m , and we assume Sk (1 < k < m) 
is the nearest sub-segment to CAR. Since Sk £ Oi and Oi ff CAR, we have Sk TT 
CAR, i.e., any point p £ s fc , it has p ^CAR. Therefore, dist(p,SP) > M PL Q, 
which means that one definitely cannot reach at the point p. Since pCsj, and 
we assume p is an any point, we can conclude that one definitely cannot reach 
at Sk- Further, since we assume Sk is the nearest sub-segment to CAR, based on 
analysis geometry, we have that one definitely cannot reach at any other sub- 
segment s r (r 7^ k and 1 < r < m), either. In summary, any sub-segment of Oi 
makes no impact on the AR. □ 



For example, line segments (h,i) and (o,p) in Figure 3(a) should be pruned 
since (M)fT CAR and (o,J?)TT CAR, too. 

Lemma 2. Given there are a number N of obstacles in 2D plane, let circle 
centering at SP with radium MPL be the CAR. If an obstacle Oi intersects 
with the CAR, and assume Oi is subdivided into a number of sub-segments by 
the crossing point(s), then, all sub-segment(s) that is/are not inside the CAR 
can be safely pruned. 

Proof. We only need to prove any sub-segment, say su CAR), makes no 
impact on AR, the proof can be derived based on lemma [1] □ 



Taking line segment (a, c) in Figure 3(a) as an example, since (a, c) intersects 



with CAR at point b and (a, b) f^CAR, so line segment (a, c) should be pruned 
partially, namely, sub-segment (a, b) will be removed. 

Lemma 3. Given there are a number N of obstacles in 2D plane, let circle cen- 
tering at SP with radium MPL be the CAR, and given an obstacle Oj divides 
the CAR into two parts, Parti and Parti, and assume SP £ Parti. Any ob- 
stacle Oi (1 < i < JV and i ^ j) can be safely pruned if it lies in the side of 
Part2. 



^ dist(p, SP) is the Euclidean distance between the two points p and SP. 
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(a) purning obsatcles (b) slicing CAR 

Fig. 3. Example of Pruning and Slicing 



Proof. Due to the limitation of obstacle Oj , one cannot enter into the side of 
Parti if (s)he do not walk out of the CAR. Otherwise, if (s)he enters into the 
side of Parti by walking out of the CAR, the travelled-distance will be larger 
than MPL, which is contrary to the given condition. Therefore, any obstacle 
that is locating in the side of Parti makes no impact on the AR. □ 



As an example, line segment (d, g) in Figure 3(a) subdivides the CAR into 
two parts, SP lies in the upper part, however, line segment (m,n) lies in the 
lower part, so (m, n) should be pruned. For the three types of obstacles discussed 
just now, we prune or partially prune according to the following method. 

First, we, using the MBR of CAR as the input, search on the constructed 
index structure, which indexing the MBRs of all obstacles, in this way, we can 
obtain the initial candidate obstacles (ICOs). As an example, the square illus- 



trates the MBR of CAR as shown in Figure 3(a) line segment (o,p) will be 
pruned since its' MBR and the MBR of CAR are disjointed, all the rest of ob- 



stacles in Figure 3(a) are ICOs. Next, we, according to their spansQ, sort all the 
ICOs in the descending order, and then deal with the sorted ICOs one by one. 
The reason we sort the ICOs is convenient for pruning correctly those obstacles 
that are relevant with the CAR but are irrelevant with the evaluation of AR. 



Taking line segment (m,n) in Figure 3(a) as an example, it lies in CAR, but it 



is irrelevant with the evaluation of AR. If we deal with (m, n) prior to (<i, g), the 
(m, n) possible cannot been pruned correctly. After we dealt with each ICOs, 
in the end, we obtain the candidate obstacles (COs). For example, the COs in 
Figure 3(a) will be (b, c), (j, k), and (e,/). 

The algorithm for pruning unrelated obstacles is illustrated in Figure Note 
that, when we deal with the sorted ICOs one by one (line 3-11), we differenti- 
ate two cases. 1) All obstacles that have been processed cannot subdivide CAR 
into two parts (line 4-5). 2) There is at least one obstacle that has been pro- 
cessed and can subdivide CAR into two parts (line 6-11). The subprocedure 



if a lino segment, say (a, fc), is parralcll to x-axis (or y-axis), the span is dist{a, b). Otherwise, the 
span is MAX(||a.a; — fr.a; || , ||a.y — where a.x is the horizontal ordinate of endpoint a. 
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Procedure Pruning.Obstacles { 
Input: J2?=l Oi, CAR 
Output: COs 



(1) ICOs<-use MBR of CAR to prune obstacle(s), switch <- 0, BSLs<- 

(2) sort all the ICOs based on their spans 

(3) for each obstacle Oi in ICOs 

(4) if {switch — 0) then //no obstacle ever divided CAR into 2 parts 

(5) Subprocedure Handle.Obstacle ( Oi, CAR,BSLs, switch) 

(6) else //obstacle ever divided CAR into two parts 

(7) for each BSL in BSLs / / BSL can divide CAR into two parts 

(8) if {Oi lies in unachievable part) then 

(9) remove Oi from ICOs, break; 

(10) if {Oi has not been removed ) then 

(11) Subprocedure Handle.Obstacle {Oi,CAR,BSLs, switch) 

(12) COs^-ICOs 

(13) return COs } 



Subprocedure Handle.Obstacle { 

Input: Oi, CAR, BSLs, switch // BSLs (Big Split Lines) 

(14) if {Oi and CAR arc disjointed ) then remove Oi from ICOs 

(15) else if (Oi lies in CAR) then // do nothing 

(16) else 1 1 Oi intersects with CAR 



(17) if (there arc two crossing point) then 

(18) BSL-<— Oi, BSLs<— BSLslJBSL 

(19) if {switch — 0) then switch — 1 

(20) Oi <— cut off the unrelated part of Oi // i.e., the part that outside the CAR } 



Fig. 4. Algorithm for Pruning Unrelated Obstacles 

Handle_Obstacle( ) is straightforward, only one thing is needed to mention, that 
is, 'BSLs' is used for storing those obstacles that have been processed and can 
subdivide CAR into two parts. 

4.3 Slicing and Placing 

Up to now, we obtained the COs, in this subsection, we discuss slicing and 
placing. On the whole, we divide the CAR into multiple parts (or sectors) based 
on the slicing lines (SLs) defined later, and place the COs into corresponding 
parts. 

Definition 3. Given a candidate obstacle COi, its' endpoints are COi.l and 
COi.2, respectively. We connect SP with COi.l / COi.2, and extend line seg- 
ment (SP,CO l .l) / (SP,CO l .2) along the direction of SP to CO l .l / CO t .2 
until it intersects with the border of CAR at a point, say CPj.l / CPi.2. If line 
segment (SP, CPj.lJ / (SP ,CPi.2) does not pass through other obstacles except 
the endpoint(s) of obstacle(s), then the line segment (SP, CPi.l) / (SP , COi.2) 
is a Slicing Line. 

Lemma 4. Given a candidate obstacle COi, it can produce two SLs at most, 
and SL at least. 

Proof. Since COi has only two endpoints, and based on Definition [21 the proof 
can be derived. □ 

For example, the dashed lines in Figure |3(b)| illustrate all the 8 SLs, line 
segment (a, b) produces two SLs, however, those small line segments produce 



S 




SL. For getting all the SLs, a straightforward solution is that, for each CO, we 
verify if it can produce SL(s), if so, we save the new produced SL(s), until we 
verified all the COs. Note that, though a SL cannot pass through other obstacles, 
but it can touch the endpoint(s) of other obstacle(s) as shown in Definition [3] 
So, we have 

Lemma 5. Given line segment (SP, CPi.l) is a SL, which is the extended line 
segment of (SP ,COi.\), where COi.l is an endpoint of candidate obstacle COi, 
if (SP, CPi.l) passes through an endpoint, COj.l, of another candidate obstacle 
COj, then the two candidate obstacles, COi o,nd COj, share a common SL. 

Proof. We only need to prove the following two claims. (1) COj can produce a 
SL by endpoint COj.l. (2) The SL produced by the endpoint COj.l is coinciding 
with (SP, CPi.l). 

Based on Dcfinition[3J we construct an extended line segment, say (SP,CPj . 1) , 
from (SP,COj.l). We first prove (SP, CPj.l) is coinciding with (SP, CPj.l). 
Since (SP, CPj.l) is extended from line segment (SP, COj.l), obviously, (SP, CPj.l) 
passes through the point COj.l. In addition, we know from the given condition 
that (SP, CPi.l) pass through the endpoint COj.l. So, we have that the (SP, 
CPi.l) and (SP, CPj.l) pass through two common points, i.e., SP and COj.l. 
Moreover, both (SP, CPi.l) and (SP, CPj.l) are the radium of CAR, Therefore, 
(SP, CPj.l) is coinciding with (SP, CPi.l). 

Next, we prove (SP, CPj.l) is a SL by contradiction. Assume (SP, CPj.l) 
is not a SL. However, we have proved that (SP, CPj.l) is coinciding with (SP, 
CPi.l), then (SP, CPi.l) is not a SL. It is contrary to the given condition. 
Hence, the lemma holds. □ 

LemmaO indicates that if a SL passes through an endpoint of other candidate 
obstacle, then this endpoint will be free from being verified. Taking candidate ob- 
stacle (a, b) in Figure [3(b)] as an example, it produces a SL, i.e., (SP, CP\), since 
(SP, CPi) passes through the endpoint c of candidate obstacle (c,d), Therefore, 
for endpoint c, we need not verify it in the future. In addition, from Lemma [SJ 
we have an immediate corollary below. 

Corollary 1. If a SL passes through a number K of endpoints of other candidate 
obstacles, then all the K endpoints will be free from being verified. □ 
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Procedure Slicing_Placing { 
Input: CAR, COs, SP 

Output: X)j'=f S " Sj.COs lj all COs has been put into corresponding sectors. 



(1) SLs<- 0, RLT <-<&// RLT means repetitive label table 

(2) for each CO G COs 

(3) for each endpoint CO.k of CO // k = 1,2 

(4) if (CO.k RLT) then 

(5) if (CO.k can produce an new SL) then 

(6) SLi ■<— the new produced SL, SLs-(— SLs (J SLi 

(7) if (SLi passes through other cndpoint(s)) then 

(8) RLT <r- RLTU (thc(sc) cndpoint(s)) 

(9) else // CO.k G RLT 

(10) remove CO.k from RLT 

(11) J2jl i "" Sj <- divide CAR into a number \\SLs\\ of sectors based on the SLs 

(12) for j = 1 to ||Sis|| 

(13) Sj.COs <— / / Sj.COs for storing COs that lie in sector Sj 

(14) EPj <— compute the enlarged polygon of Sj 

(15) for each CO G COs 

(16) for each EPj G rff f s " EPj 

(17) if (CO lies in EPj) then Sj.COs ^- Sj.COs\J CO; break ; 



(18) return Ejff S " -Sj.COs } 

Fig. 6. Algorithm for Slicing and Placing 



In the end, all the SLs are obtained, these SLs divide the CAR into multiple 
regions, i.e., sectors. For example, there are 8 sectors as shown in Figure [3(b)| 
Once the CAR is divided into multiple sectors, next, we place all the COs into 
corresponding sectors, the method to place COs is as follows. 

First, for each sector Sj, we compute its' enlarged polygon, EPj, such that (1) 
Sj C EPj, and (2) any COe EPj, it can ensures CO<E Sj, where 1 < j < \\SLs\\. 
Note that, an enlarged polygon with less edges will contribute to the follow- 
up calculations, we use an example, as shown in Figure [5j to illustrate how to 



obtain an enlarged polygon from a sector. Polygon in Figure 5(a) illustrates the 
enlarged polygon where Sj is a major-arc sector, similarly, polygons in Figure 



|5(b)| and Figure 5(c) illustrate the enlarged polygons where the two sectors are 



a semi-circle sector and a minor-arc sector, respectively. The reasons we obtain 
EPj is that to determine " COG EPj " is more simple than to determine "Co 
6 5/'. 

Next, for each CO, we verify the CO is locating in which u EPj" , it can be 
achieved by testing whether a EPj contains both of the endpoints of the CO. If 
so, we place the CO in sector Sj, which is corresponding to the enlarged polygon 
EPj. In this way, all COs will be placed in their corresponding sectors. 

Figure[6]depicts the algorithm for slicing and placing. Specifically, "line 1-11" 
is used for slicing the CAR, "line 12-17" is used for placing the COs into their 
corresponding sectors. Note that, 'RLT' is used for storing the endpoints that 
are sharing a common SL, as Lemma [5] claims. 

4.4 Evaluation of AAR 

It is easy to know, given we can obtain the accurate achievable region (AAR) in 
each sector, the AR can be obtained by merging all the AARs. In this subsection, 
we discuss how to obtain the AAR. 
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Generally speaking, we can classify the sector Sj into four cases according to 
the candidate obstacle(s) it contains. 

— Case 1: There is no CO locating in Sj. In this case, the AAR is equal to 
the sector Sj. For example, the AAR in S\ is equal to Si as shown in Figure 

m\ 

— Case 2: One CO lies in Sj, and both of the endpoints of the CO touch the 
border of CAR. In this case, the AAR is equal to a triangle. For example, 
the AAR in S5 is equal to A(SP, g, h) as shown in Figure [3 (b)| 

— Case 3: One CO lies in Sj, and only one endpoint of the CO touches the 
border of CAR. In this case, the AAR is equal to the union of two regions, 
which are a triangle and a small sector, respectively. For instance, the AAR 
in SV is the union of A(SP,e,f) and small sector <(e,m,n), as shown in 
Figure [3(b)] 

— Case 4: All the the rest of situations, which including: 

• One CO lies in Sj, and no endpoint of the CO touches the border of 
CAR. E.g., S 3 in Figure [3(b)] 

• There is no less than two COs locating in Sj. E.g., S2, S4, Sq and Sg in 
Figure [3(b)] 

Among the four cases, to evaluate the AAR for the first three cases is straight- 
forward as presented above, however, for case 4, it is a complicate situation es- 
pecially when there are a large number of COs in Sj. In the sequel, we address 
the method computing AAR for case 4. The algorithm for computing AAR is 
depicted in Figure [7J Simply to speaking, there are several steps. First, for each 
endpoint of each CO in Sj, we evaluate the shortest path (STP) from SP to the 
endpoint (line 2-4). Second, we find the visibility region, VRi.k, of the endpoint 
in enlarged polygon EPj if STP is less than the MPL (line 6). Third, we let the 
subtraction between MPL and STP be the rest of walking length (RWL), and 
use the circle centering at the endpoint with radium RWL to prune the visibility 
region, termed the pruned result as PVRi.k (line 7-10). Fourth, we combine all 
the pruned visibility regions (PVRs) and return the AAR in Sj (line 11-12). 

Procedure Computing_AAR { 

Input: Sj.COs, Sj, EPj, SP, MPL 
Output: AAR 



(1) AARf- // constructing visibility graph in this step will be more efficient. 

(2) for i = 1 to || Sj .COs || 

(3) for each endpoint CO;.fc of CO; // k = 1, 2 

(4) STP<— the shortest path from SP to this endpoint 

(5) if (STP< A/PL) then 

(6) VRi.k 4— evaluate the visibility region from the endpoint, 

(7) RWL<— A/PL-STP 

(8) TempC circle centering at COi.k with radium RWL 

(9) PVRi.k 4— VRi.k D TempC // prune those unachievable region 

(10) put PVRi.k into PVRs // PVRs is used for storing 

(11) AAR<— computing the union of all PVRs stored in PVRs 

(12) return AAR } 



Fig. 7. Algorithm for Computing AAR 



11 



Lemma 6. Given any sector Sj, which consists of a circular arc and two SLs, 
say SLi and SL2, and given we have obtained the AAR in Sj, say AARj, then, 
no matter how complicate the boundary of AARj is, it always include two edges, 
SLi and SLi- 

Proof. If the Sj is belonging to case 1, or 2, or 3 sector, it is obvious as shown 
in Figure [3(b) | Even if the Sj is belonging to case 4 sector, any point on the two 
SLs is achievable based on Definition |3l so, any point on the two SLs must be 
included in AARj. In addition, since AARj C Sj always hold, and SL\ and SL2 
are the boundary of Sj. Then, SL\ and SLi must be the boundary of AARj. 
Thus, the lemma holds. □ 

From Lemma |6l we have an immediate corollary below. 

Corollary 2. Given any two adjacent sector Sj and S* , they share a common 
SL, say SL\, the AARs in Sj and S* are AARj and AAR* , respectively. Then 
AARj and AAR* also share a common edge, this edge just is the SL\. □ 

In above algorithm (cf. Figure[7]), there are two sub-tasks, "calculating short- 
est path from SP to the endpoint COi.k" and "evaluating visibility region from 
the endpoint COi.k in EPj" , which are the key of points. In order to better 
understand the two sub-problems, we explain the details in the sequel. 

Visibility region of endpoint COi.k in EPj one classical and easy-to- 
implement method for computing visibility region is ray-sweeping. Though most 
existing algorithms for finding visibility region are focusing on the situation 
where the obstacles are polygons, their ideas can be easily incorporated into this 
sub-problem. First, we draw ray from endpoint COi.k toward the endpoint, say 
COj.k, of other line segments that are locating in EPj. For clearness, we call 
this ray as RmOi.k,CO •.&)• There are three cases as follows. 

— Case 1: The ray passes through nothing before it reaches at COj.k, but 
passes through other line segment (s) Q, before it reaches the border of EPj. 
For ease of discussion, we call the first line segment the ray passes through 
as COf. 

— Case 2: The ray passed through other line segment before it reaches at the 
endpoint COj.k. 

— Case 3: The ray passed through nothing except COj.k before it reaches the 
border of EPj . 

If the ray belongs to case 2, we discard it. For case 1, we compute the crossing 
point between the ray and COf, and store this crossing point. For case 3, we 
compute the crossing point between the ray and EPj, and store this crossing 
point. Moreover, for all 'case 1' and 'case 3' rays, we sort them in the direction 
of counter-clockwise w.r.t. polar plane coordinate. Then, the visibility region 
can be constructed based on the border of EPj and a series of crossing points 
and endpoints. We use an example to illustrate the procedure of constructing 
visibility region. 
3 

In here, if the ray only passes through the endpoint of other line segment, we consider it passes 
through nothing. That is, we do not take it as 'case 1' ray. 
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(c) (d) 
Fig. 8. Example of Visibility Region and Visibility Graph 



For example, polygon (Vi, V2, V3, V4, V5) illustrates the EPj, and the endpoint 



a illustrates the COi.k as shown in Figure 8(a) For obtaining the visibility 
region, at first, we draw ray from endpoint a toward the endpoints of other 
line segments that are locating in EPj. In here, R( a ,e) an d R(a,g) belong to the 
'case 3' ray, they intersect with EPj at V§ and V2, respectively. R( a ,c) an d R(a,d) 
belong to the 'case 1' ray, they intersect with line segment (e, /) at V7, and with 
line segment {g, h) at Vg, respectively. Given we construct the visibility region 
starting from the R( a ,e) an d traversing according to the direction of counter- 
clockwise, then, Vq first is inserted into the resultant list, next, e is inserted. For 
R(a,c)i we insert V-j followed by inserting c, since Vj and e are on the same line 
segment (e, /), as well as e has been inserted just now. Similarly, d, h, g and V2 
are inserted orderly. In this time, all rays are dealt with, we then insert the rest 
of vertexes that belong to EPj , starting from the next node of V2 , according to 
the direction of counter-clockwise, until the previous node of Vq is traversed. So 
V3, V4, V5, V\ are inserted orderly. At last, we get the visibility region, which is 
a star-shaped polygon, namely, (V^, e, V7, c, d, Vs,g, V2, V3, V4,, V5, Vi), as shown 
in Figure [8(b)] 

Shortest path from SP to endpoint COi.k in EPj Obviously, this sub- 
problem can be reduced to the shortest path problem between two points in the 
plane in which there are a lot of line segment obstacles. Lee et al. [IS] address 
shortest path among line segment obstacles but they limit to the obstacles being 
"parallel" line segments. Moreover, several papers [16117] address shortest path 
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among rectilinear obstacles in which obstacles are "rectangles" (e.g., in VLSI de- 
sign). In particular, there are a large number of literatures investigating shortest 
path among polygon obstacles in the past decades [18119120121 22 23 24]. At first 
glance, the sub-problem also is different from theirs. Though it is so, one clas- 
sical method based on graph theory can be used for tackling the sub-problem 
by minor modification |20) . On the whole, there are two main steps in order to 
tackle this sub-problem. 

— Step 1: Constructing visibility graph. 

— Step 2: To compute the shortest path based on standard shortest path 
algorithm such as Dijkstra algorithm [2"5] . 

Ray-sweeping discussed in previous paragraph also can be used for computing 
the visibility graph. For two endpoints that are from different line segments, if 
the two endpoints are visible each other, then we connect them, in this way, we 
will build the visibility graph in the end. Note that, besides the endpoints of all 
line segments that are locating in EPj, the point SP also should be constructed 
in the visibility graph. 

As an example, four line segments illustrate the obstacles, as shown in Figure 



8(c) the polygon (Vi, V<i, V3, V4, V5) illustrates EPj. Then, based on the method 
discussed above, we can get the visibility graph, in which (a,b,c,d,e,f ,g,h,V\) are 
the vertexes, line segments (a, 6), (c, d), (e, /), (<?, h) and all dotted lines are the 
edges, as shown in Figure [8(d)] 

The length of each edge in the visibility graph can be obtained by calculat- 
ing the Euclidean distance. Therefore, resorting to the Dijkstra algorithm, we 
can compute the shortest path from SP to any endpoint of these line segment 
obstacles. 

Remark 1. Since evaluating the shortest path for each pair of endpoints in Sj is 
based on the same visibility graph, we can construct the visibility graph at the 
first step in Figure [71 which will save much more time. 



4.5 Data Representation and Geometry Operations 

Data representation Up to now, we have presented the detailed algorithm for 
the FAR problem. In this subsection, we address another important issue, that 
is, how to represent the AR as well as related intermediate results such as AAR 
7 

We observe that all geometries (e.g., AR, AAR) involved in our algorithm 
can be classified into three cases in terms of their boundaries: 1) the boundary 
is consisting of only straight line segments, or 2) the boundary is consisting of 
only circular arcs, or 3) the boundary is consisting of the blend of straight line 
segment(s) and circular arc(s). Therefore, we need an unified data structure to 
represent them. In particular, this data structure should be convenient for exe- 
cuting related geometry operations (e.g., computing "AARLIPV Ri.k" as shown 
in line 10 of Figure UJ. 
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As we know, the boundaries of traditional polygon can be represented by 
a series of vertexes. However, for polygon with circular arcs (called circular- 
arc polygon), this method is invalid since two vertexes in here cannot exactly 
determine a circular arc segment. Fortunately, we observe that two vertexes in- 
corporated with an appendix point can exactly determine a circular arc segment. 
Taking the circular-arc polygon (grey region) in Figure [S] as an example, obvi- 
ously, we cannot exactly represent it if we only use 7 vertexes, V\, V%, ■ ■ ■ , V7, 
the reason is that the circular arc ViVj is ambiguity, it may be a major arc or 
a minor arc. However, when we add an appendix point Vg, then the circular arc 
ViVgVV is without ambiguity. 




Fig. 9. Example of Circular-arc Polygon 



In view of the above observation, we devise a data structure that can repre- 
sent circular-arc polygon and traditional polygon in an unified manner, in this 
paper, we use a double linked list to represent them. The most key of point, 
comparing with the representation for traditional polygon, is that we add a 'tag' 
domain, which indicates whether a point is vertex or an appendix point, other 
domains are as similar as the traditional structure, e.g., 'data' domain is used 
for storing the coordination of vertex (or appendix point) , i.e., (xi, y{), 'pointer' 
domain is used for pointing to the previous/next vertex (or appendix point). For 
clearness, in the rest of this paper, an appendix point is denoted by Vi, and a 
traditional vertex is denoted by Vi . For example, the circular-arc polygon (grey 
region) in Figure [9] can be denoted by {Vi, Vz, V3, V4, V5, V&, V7, Vg}. 

Geometry operations In the sequel, we, based on the above data representa- 
tion, address the geometry operations on circular-arc polygons. In particular, we 
mainly discuss computing the union between two circular-arc polygons, since it 
is a typical operation in the process of computing AR. On the whole, it involves 
the following two steps. 

— Step 1: Computing all the intersections between two circular-arc polygons. 

— Step 2: Constructing the resultant polygon based on these intersections and 
corresponding rules. 

A simple method for calculating the intersections is to compare all the edge 
pairs of both polygons. However, this method bears a heavily time cost especially 
when the number of edges are much more. One elegant method for reporting 
the intersections among a set of segments is the plane sweep algorithm [26127] . 
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(a) line segments (b) circular arcs 

Fig. 10. Example of Plane Sweep 



We, in here, present a modified version. The basic idea is as follows. Sorting 
the endpoints of all segments based on their x-coordinators at first. Next, the 
plane is swept (from left to right) with a vertical line. At every moment the line 
segments that intersect with the sweep-line are stored, ordered from bottom to 
top, in a data structure (e.g., a sequence linked list), DS for short. In particular, 
when the sweep line touches the left endpoint of a line segment, the line segment 
is put into DS. On the contrary, it is popped out from DS when the sweep 
line touches the right endpoint of this line segment. Therefore, the line segments 
stored in DS are changing dynamically. Note that, if the order of line segments 
in DS is changed when the sweep line moves from left to right, it must produce 
the intersection(s)|3 So the intersections among line segments can be detected 
by observing whether the order of line segments (in DS) is changed, this is the 
main difference. (Note that, original algorithm always test if a new inserted line 
segment intersects with the above/below line segment. Analogously, when a line 
segment is popped out from DS, original algorithm also do this kind of test, the 
time for both of the tests can be save using the modified version.) 

For example, on the 4 th sweep line, line segments 11 A, C, B, D" are sorted (in 



DS) based on their y-coordinators, as shown in Figure 10(a) Next, on the 5 th 
sweep line, the order of these four line segments takes place a change. Specifically, 
line segment D is promoted in two levels, so "D" must intersect with U B" and 
"C" in the interval between the 4 th and 5 th sweep line. By the way, "D" is 
popped out from DS after the 5 th sweep line, then, the sequence is 11 A, C, B" 
now. On the 6 th sweep line, the sequence is changed to "A, B, C" , we, similar 
to previous reason, can know that "£?" must intersects with "C" in the interval 
between the 5 th and 6 th sweep line. 

In here, we resort the principle of plane sweep method for finding the inter- 
sections between two circular-arc polygons. It is worth noting, since a circular 
arc possible intersects with the vertical line (or sweep line) at two points, it will 
make the plane sweep algorithm invalid. For tackling this problem, we, at first, 



^ The complete proof can be found in literature 1261 . 
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decompose non-x-monotone circular arc (if exist) into two x-monotone circular 
arcsH Specifically, we use a horizontal line passing through the center of corre- 
sponding circle to decompose a non-x-monotone circular arc into two x-monotone 
arcs. For convenience, called them as upper arc and below arc, respectively. 

Taking the big circular arc in Figure 1 10 (Id) | as an example, which is a non-x- 
monotone arc, we decompose it into two arcs 11 B" (upper arc) and "C" (below 
arc). On the 2 nd sweep line, the sequence stored in DS is "B, C, A" . For the 3 rd 
sweep line, before putting "D" into DS, the sequence has changed to "B, A, C" , 
so we know "A" must intersects with "C" in the interval between the 2 nd and 
3 rd sweep line. As the sweep line moves from left to right, similarly, we can know 
"D" intersects with "B" , and "A" intersects with "B" . 

Once the intersections between two circular-arc polygons are obtained based 
on the above method, we then construct the resultant polygon. For traditional 
polygons, one classical method constructing resultant polygon is based on the 
so called entry-exit rules |28l29j . Their basic idea is as follows. Inserting the 
intersections into the two polygons, and marking these intersections with entry 
or exit property, next, traversing from a starting vertex along the edges, in which 
it shifts to traverse another polygon once it meets the entry or exit point, and 
vice verse. In the end a circuit is produced, this is just the boundary of resultant 
polygon. 

In here, we extend this idea to our problem, but we have to do some additional 
efforts since the circular arcs present. The first effort is to sort intersections on a 
circular arc according to the direction of ant i- clockwise, this is ready for inserting 
new appendix point and for traversing. The second effort is to insert new appendix 
point into the circular arc if there is(are) intersection(s) on it, this is ready for 
eliminating the ambiguity. 

Lemma 7. Given there are a number K of intersections on a circular arc, then 
it need insert at least K and at most K + 1 new appendix points for eliminating 
the ambiguity. 

Proof. Since K intersections can subdivide a complete circular arc into K + 1 
small circular arcs, and for each small circular arc, one appendix point is needed 
and is enough to eliminate the ambiguity. So for K + 1 small circular arcs, 
K + 1 appendix points are needed and it is enough to eliminate the ambiguity. 
Further, since there is an appendix point beforehand. Therefore, when there is 
no intersection coinciding with this appendix point, only K new appendix points 
are needed. Otherwise, K + 1 new appendix points are needed. □ 

For example, there are two circular-arc polygons, as shown in Figure 1111 
CAP 1 = {Vi,V 2 ,V 3 ,Vi,V 5 ,V 6 ,V^,V 8 ,V 9 ,V w } a nd CAP 2 = {V n ,V 12 ,V^,V u ,V 15 , 
Vie, V17, Vis, Vig, V20}, in which we assume we have obtained the intersections 
between the two circular-arc polygons based on sweep line method. The black 
filled circle, black hollow circle, grey filled circle and grey filled circle with black 

a x-monotonc circular arc is such a circular arc that is at most one intersection with any vertical 
line. 
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Fig. 11. Example of Constructing Resultant Polygon 



edge denote vertexes, appendix points, intersections and new appendix points, 
respectively. For CAP\, on the arc ^V7^, there are two intersections Iq and 
I 7 , we sort them according to the direction of anti-clockwise. In particular, we 
view this circular arc as three small circular arcs which are separated by the 
two intersections. So we insert two new appendix points N\ and N 2 for eliminat- 
ing ambiguity. Similarly, for CAP 2 , we view V12V13I4 as two small circular arcs 
which are separated by the intersection I3, and we also insert a new appendix 
point N3. 

After intersections and new appendix points are orderly inserted into the two 
initial circular-arc polygons separately, we get two new linked lists, then, we as- 
sign the intersections with Entry or Exit property alternately. At last, we traverse 
the two new lined lists and obtain the resultant polygon. For instance, in Figure 
fTTl the new linked list (corresponding to CAPi) is {Vi, ^3, V2, h, V3, V4, Vg, I5, Vq, 
N2, 16, Vj, I7, Ni, Vs, V9, Ig, Vio, h, I2}, another new linked list (corresponding to 
CAP 2 ) is {Vix, l2,V 12 ,Vu,I 3 ,W 3l V u ,h,h,V lb ,h, V 16 ,I 7 ,V 17 ,Ws,V 19 ,V 20 ,I s , 
Ii}. We assign "1%, I 2 , ■ ■ ■ , Is" with "entry, exit, ■ ■ ■ , exit" , respectively. Then, 
we traverse the first new linked list, starting from V\, along the direction of 
anti-clockwise. Once we meet a intersection with "entry" property, we shift to 
traverse the second new linked list. Similarly, once we meet a intersection with 
"exit" property (in the second new linked list), we shift to traverse the first 
new linked list, and so on. In the end, we get the resultant polygon, {V\, I3, N 3 , 
V u ,h,V 3 ,%, V b ,h,V lb ,h,V 7 , I 7 ,V 17 ,V^, V 19 ,V2D,h,V 10 ,Ii,Vn,l2 }, which 
is just the union of GAP\ and CAP 2 . 
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In the process of computing achievable region, there are some other geometry 
operations, which possibly don't aim at two circular-arc polygons, or possibly 
is not computing the union of two geometries. In addition, some degenerations 
possible happen. For these issues, the complete description is beyond the scope 
of this paper, these can be either found in literatures |28I29| . or easily deduced 
based on the above discussion. We, for saving space, don't further talk the details. 

5 Correctness and Complexity 
5.1 Proof of correctness 

For proving the correctness of our proposed method, we only need to prove the 
following two propositions. 

Proposition 1. Any possible location (or point) at which one can reach is lo- 
cating in AR. 

Proof. By contradiction. Given a point is achievable, we assume it is not locating 
in AR. Obviously, this point should be either outside the CAR or inside the CAR. 
If it is the former, the point must be unachievable since the MPL is the radium 
of CAR, it is contrary to the given condition. If it is the latter, this point should 
be inside one of sectors, Sj. Recall the algorithm Framework_FAR (cf. Figure [2]), 
the AR is the union of all AARs, i.e., AR=|J^ sil AARj, where AARj is the 
AAR in Sj. Further, we know AARj C Sj. So this point only can locate in the 
region "Sj - AARj". 

Since we classify Sj into four cases (cf. Section |4~4|) . For case 1, the region 
"Sj — AARj" is equal to 0. Obviously, this point cannot locate in this kind of 
sectors. We assume this point lies in so called case 2 sector, the region "Sj — 
AARj" don't equal to 0, but any location inside the region "Sj — AARj" is 
obviously unachievable because of the limitation of MPL and the impact of 
obstacle (cf. Figure |3(b)[ ), it is contrary to the given condition. Further, we 
assume this point lies in so called case 3 sector, we also can easily and similarly 
derive that this point is unachievable, it is contrary to the given condition. 

At last, we assume this point lies in so called case 4 sector, the region "Sj — 
AARj" also don't equal to 0. Recall our algorithm computing_AAR (cf. Figure 
O, the AAR in Sj, i.e., AARj, is the union of all PVRi.k. Since PVRi.k is 
the pruned visibility region, that is, the region that the point visibility region 
VRi.k (in polygon EPj) subtracts the circle centering at the endpoint COi.k 
with radium RWL, this ensures that all locations at which one can directly 
reach from the endpoint COi.k within distance RWL are included in AARj. In 
addition, since we compute the PVRi.k for each endpoint in Sj only if RWL 
larger than 0, where RWL is the value that MPL minuses the shortest path 
from SP to COi.k, these ensure that all locations (in Sj) at which one can 
reach are included in AARj. Since we assume this point locates in the region 
"Sj — AARj" , but any location in this region is obviously unachievable, it is 
contrary to the given condition. Hence, the proposition holds. □ 
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Proposition 2. No exist a location (or point) at which one cannot reach is 
locating in AR. 

Proof. By contradiction. Given there is a possible point at which one cannot 
reach, we assume it is locating in AR. Obviously, this point should be inside the 
CAR since ARC CAR. More clearly, this point should locate in certain AARj ( 
C Si). Since we classify Sj into five cases (cf. Section I4T4")) . then, we, first of all, 
assume this point locates in so called case 1 sector (cf. Figure |3(b)[ ), obviously, 
any point in this region is achievable, it is contrary to the given condition. Next, 
we assume this point locates in so called case 2 sector, the triangle region is 
AARj, obviously, any point in the triangle is achievable, it is contrary to the 
given condition. Further, we assume this point locates in so called case 3 sector, 
we also can easily and similarly derive that any point in AARj is achievable, it 
is contrary to the given condition. 

At last, we assume this point locates in so called case 4 sector. Since PVRi.k 
is the region that the point visibility region VRi.k (in polygon EPj) subtracts 
the circle centering at the endpoint COi.k with radium RWL, and since RWL is 
the value that MPL minuses the shortest path (from starting point SP to this 
endpoint), these ensure any point in region PVRi.k is achievable. In addition, 
we know AARj is the union of all PVRi.k (cf. Figure [7]), so any point in AARj 
is achievable. Since we assumed this point locates in AARj, then, this point is 
achievable, it is contrary to the given condition. Hence, this proposition holds. □ 

5.2 Computational complexity analysis 

The algorithm framework consists of six steps (cf. Figure [2]), we analyse the com- 
putational complexity on each step, then get the total computational complexity 
by accumulating. For convenience, in the rest of the paper, we use complexity 
and computational complexity interchangeably, and we denote the complexity 
for initializing (line 1), pruning unrelated objects (line 2), slicing and placing 
(line 3), computing the AAR in Sj (line 5), computing the union of all AARs 
(line 6), and returning result by C mlUa i, C prue , C sp , C aar , C union and C return, 
respectively. So, the total complexity, Cfar, can be derived as follows. 

CfAR — Ci n m a l -\- Cprune C ' S p -\- h * C a ar C U nion -\- Creturn (^'-0 

where k is the number of sectors, it is equal to ||SXs||. Obviously, both Ci n itiai 
and C re turn are O(l) time. In addition, in the worst case, there are a number N 
of COs in CAR, so we have 

Corollary 3. Given all the N obstacles are COs, namely, all obstacles lie in 
CAR, then there are at most 2*N slicing lines. 

Proof. Based on the condition ||COs|| = N and Lemma |4j the proof can be 
derived. □ 

Therefore, the Equation [5J] can be rewrote as follows. 

CfAR — Cp rune + Csp + m«Il<K2A'{^ * C a ar} + C U nion (5-2) 
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Complexity for pruning unrelated obstacles (C prune ) Recall algorithm 
for pruning unrelated obstacles (cf. FigureE]), the subprocedure Handle_obstacles() 
takes constant time, denoted by O(ci). So, at the worst case, that is, when switch 
always is equal to (line 4), the for circulation (line 3-11) needs 0(fci *ci) time, 
where k\ = ||ICOs||. In addition, classical R-tree for finding ICOs (line 1) takes 
0(N) time at the worst case, and sort the ICOs (line 2) using classical Heapsort 
algorithm takes 0{k\ * log fci) time at the worst case, and returning COs (line 
13) takes O(l) time. Thus, we have 

Cprune = 0(N + fci * fog fci + fci * Cl) (5.3) 

In particular, ||ICOs||=N in extreme case, so Equation 15.31 can be rewrote as 
follows. 

Cprune = max < kl <N {0(N + fci * (ci + log fci)} (5.4) 

Complexity for slicing and placing (C sp ) Recall algorithm for slicing and 
placing (cf. Figure [5]), to determine if an endpoint can produce a slicing line (line 
5) takes 0(fc2) time at the worst case, where fc2 = |jCOs||. In addition, at the 
worst case, "CO.fc ^ RLT" always holds (line 4), namely, RLT always is equal 
to 0, so we have to determine (for each endpoint) if it can produce a slicing line. 
Therefore, to find out all slicing lines (line 2-10) takes 0(2 * fc 2 * fc 2 ) time, i.e., 
(9(fc2 2 ) time. Dividing the CAR into multiple sectors (line 11) only need sort the 
slicing lines w.r.t. polar plane coordinate, which takes 0{k * log fc) time, where 
fc = ||SLs||. To compute the enlarged polygon of Sj takes constant time, denoted 
by 0(c2), so "line 12-14" takes 0(k * c^) time. Moreover, since EPj at most 
have 7 edges as shown in Figure[5l to determine if CO lies in EPj (line 17) takes 
constant time, denoted by 0(03). Therefore, placing COs in Sj (line 15-17) takes 
0(fc2 * fc * c 3 ) time. Last, returning Y^j= 1^ Sj.COs takes O(l) time. Then, we 
have 

C sp =0(k2 2 + fc * log fc + fc * C2 + fc 2 * fc * C3) (5-5) 

In extreme case, ||COs|| = ||ICOs|| and ||ICOs||=N, at the same time, based on 
corollary [21 the Equation 15.51 can be rewrote as follows. 

C sp =max o<k 2 <N.o<k<2N{0(k 2 2 + k* (log A; . 

, \ . 7 7 n ( 5 - 6 ) 
+ c 2 ) + fc2 * k * c 3 )} 

Complexity for computing the AAR (C aar ) Recall algorithm for computing 
AAR (cf. Figure [7]), we construct visibility graph at the first step (as Remark 
[1] suggests) for improving efficiency, using ray-sweeping method to construct 
visibility graph takes 0((2 * fc 3 + l) 2 * log (2 * k 3 + 1)) timeQ, i.e., 0(k 3 2 * logfca) 
time, where k$ = \\Sj.COs\\. 

In addition, to compute the shortest path (line 4) using Dijkstra algorithm 
takes 0((2 * fc3 + 1) 2 ) timed, i-e., 0(k^ 2 ) time. For obtaining the visibility region 

^ This can be derived based on literature 1201 . "(2 * ^3 + 1)" is corresponding to the "n" in their 

paper, where n is the number of obstacle corners. 
^ More details about computing the shortest path please refer to the book 1301 . 
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(line 6), we first find out so called case 1 and 3 rays, it takes <3((2*(£; 3 — l)) 2 ) time. 
At the worst cases, all rays are so called case 1 and case 3 ray, sorting them takes 
(3(2*(fc 3 — l)*log(2*(fc 3 — 1))) time, next, constructing visibility region based on 
crossing points and endpoints takes 0(2 * (fc 3 — 1)) time. Therefore, to compute 
the visibility region takes 0((2*(fc 3 -f )) 2 +2*(fc 3 -f )*log(2*(fc 3 -f ))+2*(fc 3 -f )) 
time, i.e., 0(fc 3 2 + fc 3 * log fc 3 + fc 3 ) time. 

For computing "VRi.k n tempC" (line 9), we need compute the intersec- 
tion^) (between each edge of VRi.k and the circle tempC), and insert the inter- 
section^) into corresponding edge. Specifically, once we obtain intersection(s) 
between an edge and tempC ', we insert the intersection(s) into this edge, then 
to deal with next edge, in this way, all the intersections can be orderly inserted 
into VRi.k after we dealt with all the edges of VRi.k. At the worst case, the 
VRi.k has a number 2 * (fc 3 — 1) + e of edges, and each one at most intersects with 
the tempC at two points, where e is a small integer denoting the edges coming 
from EPj. So, to compute the intersection(s) and insert them into VRi.k takes 
(9(fc 3 ) time. Before we construct the result of "VRi.k O tempC" , we also have 
to sort all the intersections on the circle tempC w.r.t. polar plane coordinate, 
it takes 0(fc 3 * logfc 3 ) time, in addition, we insert the appendix point between 
each pair of intersections on the circle, it at most takes 0(fc 3 ) time. Next, to 
construct the resultant polygon based on " entry -exit" rule, it needs traversing 
all endpoints (including so called entry-exit points), it at most takes 0(fc 3 ) time. 
Therefore, to compute "VRi.k n tempC" takes 0(fc 3 + fc 3 * log fc 3 + fc 3 + fc 3 ) time, 
i.e., 0(fc 3 + fc 3 * logfc 3 ) time. 

At the worst case, "STP < MPL" always holds (line 5). Therefore, "line 
2-10" takes 0(2 * fc 3 * (fc 3 2 + (fc 3 2 + fc 3 * log fc 3 + fc 3 ) + (fc 3 + fc 3 * log fc 3 ))) time, 
i.e., 0(fc 3 3 + fc 3 2 * logfc 3 ) time. 

Next, we analyse the complexity for computing all the union of pruned vis- 
ibility region (line 11), which can be reduced to computing the union of two 
circular-arc polygons as discussed in Subsection 14.51 

Theorem 1. Given two circular-arc polygons with m and n edges, respectively, 
and assume there is no non-x-monotone arc in the two polygons, then, to compute 
the intersections using the modified plane sweep method takes 0((m+n)*\og(m+ 
n) + i) time, where i is the number of intersections. 

Proof. The proof is analogous to the argument in Theorem 3 of literature [57] , 
we omit the details due to space limit. □ 

Since a circular-arc polygon with m edges at most has a number m of non- 
x-monotone arcs, and to decompose a non-x-monotone arc takes constant time, 
then, based on Theorem [TJ we have an immediate corollary below. 

Corollary 4. Given two circular-arc polygons with m and n edges, respectively, 
then, to compute the intersections using the modified plane sweep method takes 
0((m+n) + ((m+n)*log(m+n)+i)) time, where i is the number of intersections. 

□ 
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To compute the union of all pruned visibility regions (PVRs), we use divide 
and conquer. At the worst case, U STP < MPV always holds (line 5), so, it 
runs log(2 * £3 + 1) turns. 

For each turn, since the number of edges in each PVR at most is equal to 
4* (k 3 — 1), the number of all intersections (among all edges from two PVRs) are 
(4 * (ks — l)) 2 at moslQ. Then, based on Corollary [U we know that computing 
the intersections in each turn takes 0((8 * (k 3 — I) + 8 * (£3 — 1) * log (8 * (fc 3 — 
l)) + (4*(fc 3 -l)) 2 )*(2*fc 3 + l)/2) time, i.e., 0{k 3 3 + k 3 2 * log k 3 + fc 3 2 ) time. In 
addition, to insert the number (4 * (k 3 — l)) 2 of intersections into the two PVRs 
takes 0(k 3 3 ) time. At the worst case, all the intersections on an arc, then sorting 
these intersections takes 0(k 3 2 *\ogk 3 2 ) time. Moreover, based on Lemma[71 to 
insert new appendix points takes 0(k 3 2 ) time at most. Further, since there are 
at most 2 * 4 * (fc 3 — 1) + (4 * (fc 3 — l)) 2 endpoints (including the intersections 
and vertexes), to construct the resultant polygon based on "entry-exit" rule, it 
takes 0(k 3 2 ) time. 

Therefore, to compute the union of all PVRs takes 0(log(2 * £3 + 1) * ((fc 3 3 + 
fc 3 2 * log k 3 + k 3 2 ) + k 3 3 + k 3 2 * log k 3 2 + k 3 2 + k 3 2 )) time, i.e., Q(logfc 3 * (fc 3 3 + k 3 2 
*(\ogk 3 + log A3 2 )) time. At last, returning the AAR (line 12) takes O(l) time, 
Thus, we have 

Caar =0(k 3 2 * log fc 3 + (k 3 3 + fcg 2 + fcg 2 * log k 3 ) 

3 2 2 ( 5 - 7 ) 

+ log k 3 * (fc 3 + k 3 * (log fc 3 + log k 3 ))) 

Similarly, at the extreme case, \\Sj.COs — ||COs|| and ||COs||=N, so Equation 
15.71 can be rewrote as follows. 

Caar =max 0<k 3 <N{O(kj * logfc 3 + (fc 3 3 + fc 3 2 * logfc 3 ) ^ 
+ logfc 3 * (k 3 3 + k 3 2 * (logfc 3 + \ogk 3 2 ))} 

Complexity for computing the union of all AARs (C un i on ) Since any 
two sectors at least have a common point, i.e., SP, and any two adjacent sectors 
have a common slicing line (SL). Based on this observation, we have 

Theorem 2. Given any two adjacent sector Sj and S* , and given we have ob- 
tained the AARs in Sj and S* , which are denoted by AARj and AAR*, respec- 
tively. Further, given AARj and AAR* consist of m and n vertexes ( do not 
include the appendix points), respectively. To compute the union between AARj 
and AAR* takes 0{m + n) time. 

Proof. To compute the union between AARj and AAR* is different from that 
discussed in Subsection 14.51 Based on Corollary [21 we need find out the common 
edge between AARj and AAR*, this common edge in fact is a SL. Since SP be- 
longs to any SL, and we use double linked list to store the vertexes and appendix 
points, so we need to find out the vertex "SP" from the two double linked lists, it 
takes O(m) and 0(n) time at the worst case. Then, to find out the common SL, 



Note that, the intersections do not include the vertexes from two PVRs. 
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only 4 times comparisons are needed at the worst case since each AAR has two 
SLs, it takes 4 * O(l) time. At last, we combine the two linked lists in 0(1) time 
since we only need adjust the pointer of double linked list. To sum up, to com- 
pute the union between AAR 3 and AAR* takes "0(m) + 0(ra)+4*0(l) + 0(l)" 
time, i.e., 0(m + n) time. □ 

To compute the union of all AARs, we first sort the AARs w.r.t. polar plane 
coordinate, it takes 0(k * logfc) time, where fc = ||SLs||. Then, we merge all 
the AARs using divide and conquer, it runs log k turns. For each turn, at the 
worst case, the number of vertexes for each AAR is equal to 2 * ||5j.COs|| + 3, 
then, based on Theorem [2j we know that merging in each turn takes 0(((2 * 
k 3 + 3) + (2 * fc 3 + 3)) * (k/2)) time, i.e., 0(k 3 * k) time, where k — || SLs||, and 
fc 3 = \\Sj.COsl Thus, we have 

Cunion = 0(k * log fc + fc 3 * fc * log k) (5.9) 

Further, based on Corollary 02 and \\S 3 .COs = ||COs|| and ||COs||=N at the 
extreme case, so Equation 15.91 can be rewrote as follows. 

C un ion =ma,X 0<k 3 <N,0<k<2N{O(k 3 *k* logfc)} (5.10) 

The total complexity (Cfar) For Equation 15.21 we substitute C prun e, C sp , 
Caar, C un ion with Equation 15.41 [531 15.81 and l5.101 then, we have 



Cfar = max o<fc!<Jv{0(A^ + fci * (c\ + logfci)} 

-I- max o<k 2 <N,a<k<2N{0(k 2 2 + fc * (log k + c 2 ) 

+ k 2 *k* c 3 )} + max <fc 3 <A'{ o ( fc 3 * lo S fc 3 + (h 3 (5-11) 

-I- fc 3 2 * log fc 3 ) + log k 3 * (k 3 2 * (log fc 3 + log k 3 2 ) 

+ k 3 3 )} + max o<k 3 <N,o<k<2N{0(k 3 *k* logfc)} 

By eliminating constant factor ci, ci and c 3 , the Equation l5.11l can be rewrote 
as follows. 

Cfar — max o<k<2N,o<k 1 <N,o<k 2 <N,o<k 3 <N{0(k 3 3 *logfc 3 

+ fc 3 2 * log fc 3 * (log fc 3 + log fc 3 2 ) + fc * log fc (5.12) 
+ fc2*fc + fc*fc 3 * log k + fci * log fci + N)} 

where fc = ||SLs||, fci = [[ICOs ||, fc 2 = ||COs ||, and fc 3 = \\Sj.COsW. It is worth 
noting that fc = 2N and k 3 = N cannot be satisfied simultaneously. At extreme 
case, fci = N, fc 2 = N and fc 3 = N can be satisfied simultaneously, and we 
observe that there is no any slicing line (i.e., k = 0) when fc 3 = N, at this time, 
Equation 15.121 has a maximum value, we substitute fc, fci, fc 2 , fc 3 with 0, N, N, N, 
respectively. So, we have 

Cfar =0(N 3 *\ogN + N 2 * log TV * (log N + log A^ 2 ) 

9 , (5.13) 

+ A + A * log A + N) 

Thus, by O-notation, the total complexity is 0(A 3 logiV). 
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6 Conclusion 



This paper proposes and formulates the FAR problem, we settle this problem un- 
der the assumption where the obstacles are line segments. A slice based method 
is presented, we demonstrate the correctness and the computation complexity 
of our proposed method by elaborate analysis. In future, we will investigate the 
FAR problem under the assumption where the obstacles are arbitrary polygons 
and delve into more efficient method. 



References 

1. ElGindy, H.A., Avis, D.: A linear algorithm for computing the visibility polygon 
from a point. Journal of Algorithms (JAL) 2(2) (1981) 186-197 

2. Welzl, E.: Constructing the visibility graph for n-line segments in o(n 2 ) time. 
Information Processing Letters (IPL) 20(4) (1985) 167-171 

3. Lee, D.T., Lin, A.K.: Computational complexity of art gallery problems. IEEE 
Transactions on Information Theory (TIT) 32(2) (1986) 276-282 

4. pang Chin, W., Ntafos, S.C.: Optimum watchman routes. In: ACM 
SIGACT/SIGGRAPH Symposium on Computational Geometry (SoCG). (1986) 
24-33 

5. Wang, B.: Coverage problems in sensor networks: A surveys. ACM Computing 
Surveys (CSUR) 43(4) (2011) 

6. Meguerdichian, S., Koushanfar, F., Potkonjak, M., Srivastava, M.B.: Coverage 
problems in wireless ad-hoc sensor networks. In: IEEE International Conference 
on Computer Communications (INFOCOM). (2001) 1380-1387 

7. Hwang, Y.K., Ahuja, N.: Gross motion planning - a survey. ACM Computing 
Surveys (CSUR) 24(3) (1992) 219-291 

8. Ghosh, S.K.: Visibility algorithms in the plane. Cambridge University Press, New 
York (2007) 

9. ORourke, J.: Art gallery theorems and algorithms. Oxford University Press, 
Oxford (1987) 

10. Ricci, A.: An algorithm for the removal of hidden lines in 3d scenes. The Computer 
Journal (CJ) 14(4) (1971) 375-377 

11. Eastman, J.F.: An efficient scan conversion and hidden surface removal algorithm. 
Computers and Graphics (CG) 1(2-3) (1975) 215-220 

12. Lingappan, L., Gangaram, V., Jha, N.K., Chakravarty, S.: Fast enhancement of 
validation test sets for improving the stuck-at fault coverage of rtl circuits. IEEE 
Transactions on VLSI Systems 17(5) (2009) 697-708 

13. Altshuler, D., Pollara, V.J., Cowles, C.R., Etten, W.J.V., Baldwin, J., Linton, L., 
Lander, E.S.: An snp map of the human genome generated by reduced represen- 
tation shotgun sequencing. Nature 407 (September 2000) 513516 

14. Andrews, J.H., Briand, L.C., Labiche, Y., Namin, A.S.: Using mutation analysis for 
assessing and comparing testing coverage criteria. IEEE Transactions on Software 
Engineering (TSE) 32(8) (2006) 608-624 

15. DT Lee, F.P.: Euclidean shortest paths in the presence of rectilinear barriers. 
Networks 14(3) (1984) 393-410 

16. Lee, D.T., Yang, CD., Wong, C.K.: Rectilinear paths among rectilinear obstacles. 
Discrete Applied Mathematics (DAM) 70(3) (1996) 185-215 



25 



17. Zheng, S.Q., Lim, J.S., Iyengar, S.S.: Finding obstacle-avoiding shortest paths 
using implicit connection graphs. IEEE Transaction on CAD of Integrated Circuits 
and Systems (TCAD) 15(1) (1996) 103-110 

18. Kapoor, S., Maheshwari, S.N., Mitchell, J.S.B.: An efficient algorithm for euclidean 
shortest paths among polygonal obstacles in the plane. Discrete and Computational 
Geometry (DCG) 18(4) (1997) 377-383 

19. Asano, T., Asano, T., Guibas, L.J., Hershberger, J., Imai, H.: Visibility-polygon 
search and euclidean shortest paths. In: IEEE Symposium on Foundations of 
Computer Science (FOCS). (1985) 155-164 

20. Sharir, M., Schorr, A.: On shortest paths in polyhedral spaces. In: ACM Sympo- 
sium on Theory of Computing (STOC). (1984) 144-153 

21. Mitchell, J.S.B.: Shortest paths among obstacles in the plane. In: ACM Symposium 
on Computational Geometry (SoCG). (1993) 308-317 

22. Kapoor, S., Maheshwari, S.N.: Efficient algorithms for euclidean shortest path and 
visibility problems with polygonal obstacles. In: ACM Symposium on Computa- 
tional Geometry (SoCG). (1988) 172-182 

23. Rohnert, H.: Shortest paths in the plane with convex polygonal obstacles. Infor- 
mation Processing Letters (IPL) 23(2) (1986) 71-76 

24. Storer, J. A., Reif, J.H.: Shortest paths in the plane with polygonal obstacles. 
Journal of ACM (JACM) 41(5) (1994) 982-1012 

25. Dijkstra, E.W.: A note on two problems m connexion with graphs. Numerische 
Mathematik 1 (1959) 269-271 

26. Bentley, J.L., Ottmann, T.: Algorithms for reporting and counting geometric in- 
tersections. IEEE Transaction on Computers (TC) 28(9) (1979) 643-647 

27. Shamos, M.I., Hoey, D.: Geometric intersection problems. In: IEEE Symposium 
on Foundations of Computer Science (FOCS). (1976) 208-215 

28. Greiner, G., Hormann, K.: Efficient clipping of arbitrary polygons. ACM Trans- 
actions on Graphics (TOG) 17(2) (1998) 71-83 

29. Liu, Y.K., Wang, X.Q., Bao, S.Z., Gombosi, M., Zalik, B.: An algorithm for polygon 
clipping, and for determining polygon intersections and unions. Computers and 
Geosciences (GANDC) 33(5) (2007) 589-598 

30. Cormen, T.H., Leiserson, C.E., Rivest, R.L., Stein, C: Introduction To Algorithms, 
Second Edition. The MIT Press, Cambridge (2001) 



26 



